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to  population  depletion  of  the  lowest  energy  level.  Other  effects  have 
also  been  considered.  These  are  a  finite  idler  beam  bandwidth  (simultaneous 
excitation  of  more  than  one  vibration  frequency)  a  quasi -monochromatic 
idler  beam  slightly  off  the  true  Stokes  frequency;  laser  pulse  shape; 
modification  of  the  result  for  different  rotational  levels. 

A  number  of  different  computer  codes  have  been  developed  to  solve 
the  equations  at  different  levels  of  complexity.  With  them  numerical 
investigation  of  the  effects  of  various  physical  factors  on  saturation 
was  carried  out.  Some  of  the  results  are  presented  graphically. 
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ANALYSIS  FOR  COHERENT  ANTI-STOKES  RAMAN  SPECTROSCOPY 

(CARS) 


Abstract 

Coupled  differential  equations  have  been  formulated  which  model 
the  generation  of  coherent  anti-Stokes  Raman  radiation  In  molecular 
gases.  The  equations  Include  not  only  the  basic  term  which  Involves 
the  coherent  Interaction  of  Incident  pump  and  Idler  (Stokes)  waves  to 
generate  the  anti-Stokes  radiation,  but  also  Include  additional  terms 
which  model  other  radiation  and  absorption  processes,  both  stimulated 
and  coherent,  which  occur  simultaneously.  These  are  important  in 
determining  saturation  effects  due  to  population  depletion  of  the 
lowest  energy  level.  Other  effects  have  also  been  considered.  These 
are  a  finite  idler  beam  bandwidth  (simultaneous  excitation  of  more 
than  one  vibration  frequency)  a  quasi -monochromatic  idler  beam 
slightly  off  the  true  Stokes  frequency;  laser  pulse  shape;  modification 
of  the  result  for  different  rotational  levels. 

A  number  of  different  computer  codes  have  been  developed  to 
solve  the  equations  at  different  levels  of  complexity.  With  them 
numerical  investigation  of  the  effects  of  various  physical  factors  on 
saturation  was  carried  out.  Some  of  the  results  are  presented 
graphically. 
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I.  INTRODUCTION 


Coherent  anti -Stokes  Raman  Spectroscopy  Is  a  form  of  Raman 
spectroscopy  which  has  developed  with  the  advent  of  laser  light.  In 
It  the  medium— usually  gas  or  liquid— Is  Irradiated  by  two  waves*  the 
pump  wave  at  frequency  cop  and  a  lower  frequency  wave,  the  Idler  wave, 
at  o»j.  When  their  frequency  difference  matches  a  vibrational  frequency 
b>y  characteristic  of  the  molecule  the  Idler  wave  is  then  at  the  Stokes 
frequency  <i)s  *  u>p  -  «jy.  One  result  Is  coherent  generation  of  a  new 
wave  at  the  anti -Stokes  frequency  u»a  *  2u>p  -  o>s  *  tup  +  d»y  provided 
the  Incident  wave  Intensities  are  adequate  to  excite  the  third  order 
nonlinear  terms  In  the  electric  susceptibility  of  the  gas.  For  col- 
linear  phase-matched  beams  the  Intensity  or  power  per  unit  area,  I  ,  In 
the  new  wave  Is  roughly  proportional  to  the  square  of  the  Incident  pump 
Intensity,  Ip,  times  the  Incident  Stokes  Intensity,  I$,  and  also  is 
proportional  to  the  square  of  the  distance  over  which  the  phase  matched 
Interaction  occurs.  Thus  a  very  Intense  anti-Stokes  output  can  be  gen¬ 
erated,  orders  of  magnitude  greater  than  spontaneous  Raman  scattering. 
Futhermore,  unlike  spontaneous  Raman  scattering,  it  is  directional  so 
that  all  the  anti-Stokes  radiation  can  be  readily  collected.  Because 
of  these  two  properties  as  well  as  others  the  generation  of  CARS  spectra 
offers  a  potentially  very  practical  tool  for  the  measurement  of  molecular 
species  concentrations  and  temperature  In  gases.  It  Is  especially  useful 
as  a  probe  of  combustion  gases  since  these  commonly  contain  many  par¬ 
ticulates  which  can  so  strongly  scatter  Incident  radiation  and  also, 
due  to  the  laser  heating,  radiate  themselves,  they  can  cause  relatively 
weak  spontaneous  Raman  radiation  to  be  lost  in  the  particle  generated  "noise" 


while  the  far  stronger  CARS  scattering  would  not  be.  Even  with  no  par¬ 
ticles  In  the  gas  ordinary  spontaneous  Raman  scattering  Is  In  fact  much 
weaker  than  the  Rayleigh  scattering  which  arises  from  the  gas  molecules 
and  there  Is  also  Interference  due  to  fluorescent  radiation.  Because 
of  these  considerations  AFAPL  Is  actively  laying  the  groundwork  for  the 
use  of  CARS  as  an  eventually  routine  diagnostic  tool  for  combustion 
systems.1 »2  * 3 

The  relation  between  measured  CARS  spectra  for  gases  and  the 
usually  desired  quantities;  molecular  species  concentrations  and  temper¬ 
ature  must  be  made  on  the  basis  of  theory  taking  Into  account  both  the 
particular  experimental  setup  as  well  as  modification  of  the  Internal 
molecular  structure  of  the  constituent  molecules  by  the  applied  and 
Internally  generated  electromagnetic  fields.  The  theory  also  requires  a 
combination  of  quantum  mechanical  considerations  regarding  the  populations 
of  the  molecular  vibrational  and  rotational  energy  levels  and  transitions 
between  them  and  classical  wave  propagation.  The  latter  accounts  for 
the  coherency  in  the  Interactions  between  the  various  waves  which  Is  a 
fundamental  part  of  the  CARS  generation  of  the  anti-Stokes  wave. 

By  limiting  the  theory  to  the  most  essential  physical  factors 
which  explain  CARS  a  convenient  basic  analytical  formula  for  the  CARS 
Intensity  results  which  gives  the  proportionalities  discussed  above. 

While  the  formula  has  clearly  been  of  great  value  in  a  general  under¬ 
standing  of  CARS  It  Is  Inadequate  for  the  accurate  quantitative  pre¬ 
diction  of  the  CARS  spectra  under  certain  conditions.  In  particular, 

It  Is  quite  Inadequate  to  predict  saturation  effects  which  occur  as 


Ip,  I$  or  both  are  Increased.  It  is  also  inadequate  to  treat  the 
case  when  the  idler  beam  is  not  quasi -monochromatic  but  contributes 
to  excitation  of  more  than  one  vibrational  line  simultaneously.  For 
saturation  effects  to  be  accurately  predicted  requires  the  inclusion 
of  stimulated  generation  and  absorption  of  the  pump.  Idler  and 
anti -Stokes  waves  as  well  as  stimulated  and  coherent  generation  of  the 
second  order  Stokes  wave  at  the  frequency  ux  *  “p  ■  2“v  =  2ws  -  wp 
"symmetrical"  with  the  anti-Stokes  frequency  2u»p  -  w$.  This  is  the 
so-called  "scissors"  CSRS  process.  All  these  processes  occur 
simultaneously  with  CARS. 

In  this  report  we  outline  the  derivation  of  the  governing 
differential  equations  very  concisely  since  much  of  it  directly 
parallels  results  already  in  the  literature  except  that  here  some 
additonal  physical  processes  and  refinements  are  included.  Aspects  of 
the  numerical  integration  problem  are  discussed  and  numerical  results 
presented  with  emphasis  on  effects  occurring  due  to  depletion  of  the 
population  of  the  lowest  vibrational  energy  level. 


II.  DERIVATION  OF  EQUATIONS  TO  MODEL  CARS 


A.  Physical  Basis  for  CARS 

To  describe  the  generation  of  CARS  spectra  we  consider  the 
vibrational  energy  level  diagram.  Figure  1,  showing  the  potential  energy 


Fi gure  1 .  CARS 

function  and  the  ground  state  and  two  higher  vibrational  energy  levels 
for  a  diatomic  molecule.  When  a  gas  of  such  molecules  is  pumped  by 
intense  laser  radiation  of  circular  frequency  o)p,  the  set  of 
simultaneous  coherent  photon-molecule  interactions  illustrated  in 
Fig.  1  can  be  set  up  by  coherent  interaction  of  the  electromagnetic 
and  vibrational  waves;  absorption  of  two  pump  photons,  emission  of  a 
photon  at  “s  3  “p  -  “i ,  the  Stokes  frequency  and  at  wa  =  2u>p  -  ws  *  wp  , 
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the  anti-Stokes  frequency.  If  the  gas  Is  also  Irradiated  by  an  external 

source  at  the  Stokes  frequency  the  emission  of  anti-Stokes  radiation 

Is  greatly  Increased,  primarily  because  the  nonlinear  term  proportional 
* 

to  I E  j 2  E  In  the  Induced  polarization  of  the  gas  is  enhanced.  In  the 

r  — 

CARS  experimental  systems  at  AFAPL,  the  gas  Is  irradiated  by  coll  inear 

pump  and  Idler  laser  beams  of  the  same  polarization.  These  are  focussed 

Into  a  small  region  In  the  gas  so  as  to  obtain  a  high  resolution  as  well 

as  a  high  Intensity.  In  the  focal  region  the  beams  "neck  down"  to  a  short 

cylindrical  region  In  which  they  may  be  considered  to  have  plane  phase 

fronts.  These  conditions  will  be  assumed  in  what  follows. 

Along  with  CARS  there  are  other  coherent  Raman  processes  which 

take  place;  inverse  CARS  and  coherent  secondary  Stokes  Raman  (leading 

to  CSSRS  spectroscopy)  and  its  Inverse.  These  are  illustrated  by  the 

diagram  In  Fig.  2  where  It  Is  seen  that  the  CSSRS  frequency  Ux  = 

2W  -  w  w  -  2m  .  There  are  also  stimulated  emissions  which  occur 
s  p  p  v 

Independently  of  one  another  nad  do  not  combine  coherently  as  in  CARS. 
These  are  shown  In  Fig.  3. 

These  processes  are  all  Induced  by  the  lowest  order  electric 
field  dependent— hence  nonlinear— term  in  the  susceptibility  of  the 
gas.  In  this  point  of  view  the  induced  molecular  polarization  may  be 
written  In  terms  of  the  electric  field  susceptibility  tensors  and 
x(3)  as 

".<*>  ■  x£>MO  *  WY  <' 


or  as 

P(t)  ■  P*1*  (t)  +  P(3*  (t) 


j 


(2) 


Inverse  CARS 


CSSRS 


Note:  the  populations  of  the  vibrational  level  (1)  and  ground  level  (0) 
do  not  change  because  of  these  processes. 


Figure  2.  Coherent  Raman  Processes. 
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Note:  the  populations  of  the  vibrational  level  and  ground  levels  change  as  each 
process  proceeds. 


Figure  3.  Stimulated  Raman  Scattering. 
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The  vector  £  is  made  of  the  sum  of  pump,  idler,  secondary 
Stokes  and  anti -Stokes  components  all  with  the  same  direction  of 
propagation  and  the  same  linear  wave  polarization. 

The  stimulated  processes  alter  the  populations  of  the  ground 
and  energy  levels.  Furthermore,  while  all  the  processes  shown  in 
Figs.  1  through  3  involve  only  levels  0  and  1  the  effect  on  the 
population  of  the  additional  process  shown  in  Fig.  4. 
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Figure  4.  Stimulated  Raman  Scattering  Involving  Third  Vibrational  Level. 

This  difference  Stokes  frequency,  w$2  =  wp  -  (u>2  -  Wj)  differs 
from  the  Stokes  frequency  u>s  =  wp  "  “l  by  a  small  percentage  only. 

The  difference  is  due  to  the  lack  of  symmetry  of  the  potential  function 
as  explained  in  Herzberg.4  Hence  both  processes  will  be 
stimulated  by  an  incoming  idler  beam  of  even  small  bandwidth.  For  the 
equations  to  be  presented  in  the  following  section  we  will  ignore  the 
process  in  Fig.  4.  It  will  be  considered  again  in  Section  o. 
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B.  Two-Level  Equations 

The  incident  waves  are  assumed  to  be  quasi -monochromatic 
coll inear  waves  with  the  same  wave  polarization.  Then  the  incident 
CARS  and  stimulated  waves  will  also  be  collinear  and  with  the  same 
polarization  as  the  incident  waves.  Accordingly,  the  total  electric 
field  may  be  written  as  a  scalar  sum 


A  (x.t) 
a 


exp(-iw  t  + 
01 


ik  x) 
a 


+  cc  , 

„ 


(3) 


where  ax  stands  for  pump  (p),  a2  for  idler  (I),  o3  for  anti-Stokes  (a), 
a4  for  secondary  Stokes  (X),  and  cc  means  complex  conjugate.  If  the 
idler  frequency  is  not  the  true  Stokes  frequency  “j  ^  us  =  “p  "  « 
then  by  u>  and  in  the  following  we  will  mean,  respectively, 

a  X 

2ai p  -  (Dj  and  2^  -  up.  This  is  done  to  avoid  additional  notation 
such  as  possibly  u)aI  and  u>XI.  We  are  including  the  possibility  of 
(1)^  f  a>s  to  be  able  to  evaluate  the  rate  of  drop  off  in  the  anti-Stokes 
field  with  small  deviations  from  the  ideal  resonant  condition  <dj  =  o>s. 

Substitution  of  Eq.  (2)  into  Maxwell's  equations  and  elimination 
of  the  magnetic  field  in  the  usual  way  leads  to  the  wave  equations 


lii .  oi  &  *  4s.  li-p(3)(x>t) 

3X2  C  3t2  C2  3t2 


(4) 


When  Eqs.  (3)  and  (1)  are  substituted  in  (4)  and  terms  with  a  comnon 
frequency  equated  one  finds  a  set  of  equations  for  the  Ea 


■  n  :  ;■  r,-,P  iMHU 
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32E 
_ _ a 

3X2 


+ 


n2(ri2  E 
c2  a 


p(*) 


(X) 


(5) 


These  equations  for  the  Ea  are  coupled  through  the  P^(x);  these  are 
terms  which  require  expressions  for  the  of  Eq.  (1)  before  they 

can  be  expressed  explicitly.  For  proper  quantum  mechanical  derivations 
of  the  third  order  susceptibilities  we  refer  to  Druet  and  Taran5 
where  it  is  shown  that,  for  the  waves  with  frequencies  up  - 
they  take  the  form  for  stimulated  scattering 


x 


(3)  =  N 

pa  h 


)c“ 
— ( 


da 

da 


)  - 

V1  (ay,  -  ay) 


1 


'+  ir 


v"v' 


(6) 


or  its  conjugate  for  coherent  scattering.  Here  a'  represents  Stokes, 
anti-Stokes  or  secondary  Stokes  and  the  transition  in  question  is  between 
vibrational  levels  v'  and  v"  =  v'  +1.  Nv„  and  Nyi  are  their  respective 
population  densities  and  N  =  Nv„  +  Nyt.  (da/dn)  is  the  spontaneous 
Raman  scattering  cross  section.* 


The  cross  section  used  here  is  ws/u>p  times  the  cross  section  in  Druet 
and  Taran.5  J 
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In  addition  to  the  above  considerations  we  use 


32E  32A  exp(-1(ut  +  k  x)  32E  3E 

—a  =  - 2_  .  - “  +  2ik  -  k2E 

2  3X2  av2  3X  a 


(7) 


3X 


3X 


The  explicit  equations  which  result  from  the  substitution  of  the 

forms  like  Eq.  (6)  into  (1)  and  (5)  and  then  equating  coefficients  of 
-iu>  t 

e  a  are 


3A 

— —  =  p 

3X  a 


(8) 


3F 
_ a 

3X 


-4ttv  i (F  -  if  ) 
a  a  a 


(9) 


where  a  represents  pump,  idler,  anti -Stokes  and  secondary  Stokes, 
i.e.,  p,  I,  a  and  X,  respectively,  and: 


*  inverse  wavelength  (cm) 


fa  =  12^IA^eiXAk+xap^l%i2+xIXWpeiX^ 


fx  =  12lr2vxCxpIApA!  eiXDk+4<XIpApAIAI  eiX6k  +  XixAxiAI|2] 


f!  *  12tt2Vi[x*pa*A2  e1xAk  +  (xpjl Ap|2Aj  +  xIX  1  AXl2Aj3 


(10 

(11 

(12 


fp  *  ’2*avpHxp,IM2  +  xapl%l!)Ap  *  (xp*.  *  *,p>»M 
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Here 


Ak  =  2k„  -  kT  -  k, 

p  1  a 


6k 


k  +  k.  -  kv 
p  I  X 


Ok 


04) 


Clearly,  this  set  of  equations  could  be  cut  almost  in  half  and 

hence  be  simpler  and  much  cheaper  to  solve  numerically  if  the  second 

derivatives  in  A  ,  the  aF  /ax  terms  governed  by  Eq.  (9)  were  negligible, 
o  a 

Furthermore,  boundary  conditions  on  the  F^  are  not  normally  available 

from  the  experimental  data  so  to  integrate  the  present  set  one  has  to 

estimate  the  boundary  conditions  by  an  analytic  iterative  procedure. 

This  is  given  in  the  Appendix. 

From  physical  considerations  one  expects  that  the  terms 

aF  /ax  will  have  almost  no  influence  on  the  results  and  we  have  verified 
a 

this  in  a  few  numerical  runs.  In  all  other  work  we  assumed  that 

aF  /ax  *  0. 
a 

For  efficient  computing  separate  codes  based  on  setting  (9) 

to  zero 

F  =  if  (15) 

a  a 

were  used  when  the  a2A  /ax2  terms  were  to  be  neglected. 

at 
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In  these  equations  we  must  also  bear  In  mind  that  the  population 
difference  (NQ  -  N^)  is  variable  in  time  because  of  the  stimulated 
Raman  Induced  population  changes  Indicated  In  Figure  2.  In  addition 
there  is  decay  from  any  value  Nj  >  N1Q  to  the  value  N1Q  where  N1Q  is 
the  thermal  equilibrium  value.  This  is  due  to  collisional  interactions 
in  which  the  vibrational  energy  of  the  molecules  is  converted  to 
vibrational  energy.  This  takes  place  with  a  relaxation  time  which  Is 
characteristic  of  the  molecule  In  question  as  well  as  the  other  molecules 
in  the  gas,  the  pressure  and  temperatures. 

To  derive  a  differential  equation  for  the  time  variation  of 
Ng  -  N-j  we  consider  that  the  photon  density  n  In  the  incoming  idler 
(Stokes)  beam  changes  to  n  +  (an/ax)dx  in  traversing  a  gas  sample  of 
length  dx.  The  intensity  of  the  Incoming  beam;  i.e.,  the  time 
averaged  Pointing  flux,  is 

Ij  3  lrAsAs  =  n  (16) 


so  that 


an 

ax 


aA, 


aA, 


(A  _ i  +  A  _ 1  \ 

T  VMi  ax  Mi  ax  ' 


I  ax 


(17) 


Me  have 


a  A| 

W  S  1fl[sRS 


(18) 
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Mhere  f  j  Is  given  by  the  stimulated  Raman  terms  in  Eq.  (12); 

f 1 1  $ r  s  1 2ir2  Vj  (xpj  Aj|Ap|2  +  XIXAj|Axl2}  •  (19) 

Then  we  find  that 


34  +  4  +  ,24irU  I  2 
at  T  hN  1  I1 


xIXllAxl2  * 


(20) 


and  aq  is  the  value  of  a  at  thermal  equilibrium,  r  is  the  relaxation 
time  for  transition  from  levels  v  =  1  to  v  =  0.  The  subscript  1 
indicates  imaginary  part.  Here  a  is  the  quantity 


A  = 


"l  -  Hi  -  Nq 

♦  Nq 


(21) 


To  find  Nv  we  use  the  Maxwellian  distribution 


-Ev/kT  ^ 
Nv  =  Ne  v  /  2. 


-Ev/kT 


(22) 


where 

Ey  =  hc[we(v  +  1/2)  -  weXg  (v  +  1/2)  ]  (23) 

The  quantities  a>e  and  u)exe  are  available  in  Herzber^  for  all  diatomic 
molecules. 


! 


! 
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C.  Simple  CARS  Formula 

The  dominant  term  in  the  generation  of  the  anti-Stokes  wave 
is  the  first  term  on  the  right-hand  side  of  Eq.  (10).  This  represents 
the  coherent  (CARS)  process.  To  the  extent  that  variation  of  Aj  and 
Ap  and  NvM-Nv<  (which  enters  into  Xpi  as  shown  in  Eq.  (6))  and  the 
other  terms  on  the  right-hand  side  in  the  equation  can  be  ignored  one 
then  has,  on  integrating  over  0  <  x  <  L 


Using 


Aa  ’  -12"2  VpfVp  « 


i  Ak 


sin  j  Ak 
jj-  Ak 


y- A  A 

dir  an 


(24) 


we  get 


I.  *  (  I5-  )2  02*2VpI>2  '■> 


2 

I‘P 


sin  Ak  \2 


Ak 


(25) 


One  purpose  of  this  Investigation  is  to  check  how  good  a  result  (25) 
gives  as  one  enters  parameter  ranges  where  the  assumptions  are  not 
really  valid. 


D.  Three-Level  Equations 

In  this  case  v  «  0,  1  and  Z.  Vibrational  transitions  directly 
between  levels  0  and  2  are  forbidden;  their  probability  is  so  low 
we  need  only  consider  the  01  and  12  transitions.  Then  Introducing  the 
two  idler  waves  1^  and  1^  of  slightly  differing  frequencies  equal  to 
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o r  near  to  w.  and  «_  where  a>  ®  w..  -  u>, ,  u>  =  u  -  (w  -  u>  ) 

c  cn  c«  n  *  m  Ca  n  a  .  ' 


Si  P 


s2  P  '  2  1' 


we  are  led  to  two  anti -Stokes  waves  u.  -  2u  -  u.  and  u  »  2u>  -  w.  . 

aj  p  I2  a2  p  I2 

As  In  the  earlier  work  the  u  .  are  not  the  "true"  quantities  u  +  u> 

al  Pi 

and  u_  +  (to  -  ai  )  unless  w.  =  u)T  ,  to,  =  u>T  ,  We  consider  the  two 
P  2  1  Si  Ij  s2  I2 

Incident  fields  and  Ej2  although  in  physical  practice  there  will 
likely  be  only  one  idler  beam.  The  equations  as  written  below  can  model 
this  situation  whether  the  beam  is  purely  monochromatic  or  has  a  band¬ 
width  sufficient  to  cover  <*>  and  w  . 

Si  s2 

The  equations  for  the  f  are  now  found  by  replacing  all  the  a 
subscripts  by  a.,  j  *  1  or  2  except  for  subscript  a  *  p.  The  A k,  6k 
and  Dk  must  likewise  be  subscripted  1  and  2.  Thus,  for  example, 


fT  ■  12u2vf 
1 


*  * 

X-  p  A 

j  L  j  j 


ixAk 

An2  e 


+  xpi  I Ap| 2  Aj  +  xj  .x.  I AX I  2  AI 

J  j  J  J 


and  similarly  for  f  and  f  ,  On  the  other  hand  the  Eq.  (13)  for 
aj  XJ 


fp  Is  replaced  by  just  one  equation, 


fp  =  1 2ir2vp  V  f  .  ,A  j 

j=l  -  Plj 


+  YaPj  I A3jI2)  AP  +  (xpl.  +  xap.) 

*  -IxAk, 

Aj  Aa  Ap  e  J 
Xj  aj  P 


*  ,  ^Dk 

(  A*  Ar  e  J 

'J  li 
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We  now  haye  to  obtain,  for  use  In  the 


both 


A 


(24) 


By  a  derivation  directly  analogous  to  that  used  for  Eq.  (16)  but 
involving  the  stimulated  Raman  processes  of  both  figures  3  and  4  we  can 

obtain  the  appropriate  equations.  In  terms  of  * 


A  =  A  and  A  *  (N  -  N  )/N 
1  2  2  0 


(25) 


!27t 

h~~N 


|Ap|2  (2|ASJ*  x 


.  +  [A  I2 
ii  1  s2' 


X2i} 


3  T  T 
20  21 


T 

10 


(  A  -  A  ) 
2  20 


(26) 


12* 


—  lEplM  |ASll2*lt  *  |A.2P  ) 


*  —  (  r~  >  <  -  a>»  > 

0  20  21  10 


(27) 


+ 


1 

Tio 


4 

T20 


)  (Az  -  A2o) 

t21 


♦ 

The  equations  In  terms  of  Aj  and  a2  are  much  better  suited  for  numerical 
Integration  than  equations  in  terms  of  A  and  aJ2  or  in  terms  of  Nx  and  N2. 
This  Is  because  the  rates  of  change  of  Ai  and  a2  are  relatively  slow  while 
this  Is  not  true  for  N2  and  aj2  near  t*0. 
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D.  Rotational  Levels 

The  rotational  energy  levels.  Eu  associated  with  each  vibration 
level  v  have  so  far  been  ignored.  If  we  limit  ourselves  to  the  Q  branch 
transitions  between  vibrational  states  there  is  no  change  in  the  J  value. 
For  this  case  the  previous  equations  are  readily  generalized  so  that 
they  apply  to  a  given  rotational  level  J.  Then  one  replaces  formula 
(18)  by  the  equilibrium  values  appropriate  to  the  Jth  levels: 


N(2J  +  l)gn  e 


-(Ev+E0)/(kT) 


v  ,Jo 


YTV" 

V 


GnkT/(21tBfic) 


(28) 


where  Ey  is  given  still  by  (19), 'll  =  h/(2n). 


Ej  *  Bfic  J(J  +  l)/kT  (29) 

and  g„  and  G„  are  constants  characteristic  of  the  molecules  in  question. 

Jn  n 

(see  Herzberg6). 

E.  Measured  Quantity-Energy 

The  main  physical  quantity  which  is  measured  experimentally  is  the 
anti -Stokes  energy  or  power  per  unit  area  incident  on  a  receptor  surface 
at  the  end  of  the  sample,  integrated  over  the  pulse  time  and  the  detector 
surface  area.  We  are  assuming  uniform  plane  waves  so  the  area  integration 
is  trivial.  Hence  the  quantity  we  shall  compute  to  model  the  experiment 
will  be  energy /area 


III.  COMPUTATIONS 


A.  Computer  Codes 

The  equation  sets  of  the  previous  section  were  programmed  in 
FORTRAN  IV  for  computation.  For  simplicity  and  speed  in  the  computations 
separate  programs,  CARSO  through  CARS5,  were  written  for  different  cases. 


Code 

Name 

Main 

Equations 

Number  of 

Energy  Levels 
(including  qround) 

Is  CSSRS 
Included? 

Are  second 
derivative 
terms  included? 

Integration 

Method 

CARSO 

(15),  (20) 

2 

no 

no 

RK 

CARS2 

(22)-(2 7) 

3 

no 

no 

RK 

CARS3 

(15),  (20) 

2 

no 

no 

A,M  and  G 

CARS4 

(8),  (9), 

(20)  2 

no 

yes 

RK 

CARS5 

(15),  (20) 

2 

yes 

no 

RK 

Here  RK  stands  for  4th  order  Runge-Kutta,  AM  for  Adams-Moulton  and  G  for 
Gear's  method.  Provision  was  made  in  CARSO  to  set  individual  terms  equal 
to  zero  on  the  right.  Also  CARS5  was  set  up  to  compute  either  Gaussian 
or  rectangular  shaped  pump  and  idler  pulses  whereas  the  other  codes  are 
limited  to  rectangular. 

Numerical  comparisons  between  sample  CARSO  and  CARS3  results 
showed  no  reason  to  favor  AM  or  G  over  the  RK  method  so  RK  was  used  for 
all  the  data  presented. 
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The  Integration  in  space  and  time  was  organized  to  proceed  between 


the  discrete  point  mn  shown  in  Figure  5  in  a  manner  designed  to 
minimize  storage. 


//  /2  /J 

'  -4.  i 


2!  22  2 3  • 

&  32  33  • 


•  *  *  • /N 

■4— * - 1— l - 


M/\  M2  M3 


‘2/V 

•J/V 


MN 


Fig.  5  Diagram  to  explain  computation  scheme. 


x 


To  describe  it  let  the  symbolic  statement  nm  -*•  (n  +  l)m  mean  to  carry 
out  the  integrations  of  the  Aq  with  respect  to  x  between  the  points  nm 
and  (n  +  l)m  while  nm  -*■  n(m  +  1)  means  to  integrate  the  equations  for 
A(or  Aj  and  a2)  with  respect  to  time  between  nm  and  n(m  +1).  The 
integrations  between  the  points  were  carried  out  using  many  further 
subdivisions  whose  size  is  adjustable  automatically  to  fit  prestated 
accuracy  requirements.  Distance  x  ranged  from  0  to  L  while  time  was 
measured  from  0  to  T,  with  L  and  T  chosen  to  correspond  to  values  ty¬ 
pical  of  experiments  at  AFAPL. 


The  computation  was  carried  out  using  the  following  order  and 
approximations: 

11  -►  21  using  a0  as  initial  value 

11  12  using  a0  in  computing  x's 

12  -*•  22  using  a0  as  initial  value 


l(n  -  1)  IN  using  a0  in  computing  x's 

Then 

21  -*•  31  using  a (21 )  as  initial  value 

21  ■*  22  using  linear  interpolation  between  (21)  and  (22)  in  computing  x's 

22  32  using  a(22)  as  initial  value 


(M  -  1 )N  MN  using  A((m  -  1 )N)  as  initial  value. 

Numerical  comparison  of  simple  CARSO  and  CARS4  results  showed, 
as  expected,  almost  no  difference  in  the  numerical  results  for  the  elec¬ 
tric  fields  Ea  or  energies  £  as  computed  with  the  second  derivative  terms 
or  neglecting  them.  Keeping  them  however  led  to  so  tremendous  an  increase 
in  computing  time,  that  it  was  Impractical  from  a  cost  standpoint  to  ex¬ 
tend  the  computations  as  far  in  distance  and  time  as  would  correspond  to  the 
CARS  experimental  work  at  WPAFB.  In  light  of  these  results  all  subsequent 
work  was  done  without  second  derivative  terms. 

Another  comparison  which  was  made  was  to  run  the  ame  data  using 
CARSO  and  CARS5.  The  results  showed  that  there  was  practically  no  change 
1n£„;  when  CSSRS  energy  developed  it  was  not  at  the  expense  of 

a  X  a 


B.  Presentation  of  the  Results 

£.  for  the  various  cases  to  be  described  has  been  plotted  vs. 

a 

I Apl **l Axl 2 ^  (2n/c)3  Iplj^on  log  paper. 

If  the  simple  CARS  formula  (25)  were  to  hold  then  the  computed 
curves  would  simply  be  straight  lines  at  45°  to  the  axes.  These  lines 
would  intercept  the  vertical  axis  at  a  height 

sin  (*-Ak)  2 

log  (f^r  )  (12n2  vaxpI)2  (L/^)Ak'  L2T  . 

As  I2It  increases  the  effect  of  the  non  CARS  terms  becomes 
non-negligible  and  with  that  goes  the  fact  that  changes  in  the  energy 
population  levels  become  important.  The  Ta  vs.  |Ep|4|Ej|2  plots  bend 
below  the  straight  line  which  results  when  these  phenomena  are  neglected. 
This  bending  is  the  saturation  effect  referred  to  in  the  abstract  and 
introduction. 

To  show  how  the  bending  of  the  S.  curves  is  tied  in  with  energy 

Ql 

level  depletion  the  corresponding  curves  of  a  or  Aj  and  a^  vs  ( Ep ( ( E j 1 2 
are  also  plotted. 

C.  Data 

Computations  were  made  for  N2.  The  corresponding  constants 
needed  to  compute  aq  and  the  x  were  taken  from  Hertzberg4.  Relaxation 
times  were  taken  from  White  and  Millikan7.  In  CARS2  we  need  several 
relaxation  times.  We  could  find  no  such  detailed  data  and  have  simply 
used  a  common  t  value  equal  to  the  above.  The  actual  data  used 


for  the  computational  results  presented  in  Section  E  is,  (in  Gaussian 


units) 

0.2331  E04 

Notation: 

Vv  = 

(E04  =  104,  etc.) 

VI  = 

0.1207  E05 

VP 

0.1440  EOS 

da/dfl  = 

0.3300  E-31 

Case  A: 

P  = 

0.7600  E  +  05 

T  = 

0.3000  E  +  04 

T  = 

0.400  E  -  06 

(Ao  ' 

■  0.5080  E  +  00) 

Case  B: 

P  ■ 

0.7600  E  +  05 

T  * 

0.1000  E  +  04 

T  = 

0.1000  E  -  04 

' 

■  0.9328  E  +  00) 

These  p,  T 

and  t  combinations 

are  from  the  White  and  Millikan^  empirical 

data  and  were  chosen  to  illustrate  the  effects  of  varying  t.  The  aq  values 
however  also  vary  and  their  computed  values  are  shown  in  parentheses. 

D.  Discussion  of  Results 

Aside  from  the  results  discussed  in  the  previous  sections  which 
do  not  lend  themselves  to  graphical  display  we  present  a  series  of 
curves  where  I^Ij  is  varied  by  changing  Ip  while  keeping  Ij  fixed. 

Figure  1  illustrates  the  buildup  of£_  with  time  by  showing 

Q 

curves  for  T  *  10  ns  and  T  *  20  ns  for  each  fixed  level  of  pump  and  idler 
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Irradiation.  The  curves  show  that  once  saturation  due  to  increasing 
Ip  becomes  strong  very  little  or  no  additional  anti-Stokes  power  is 
delivered  in  the  last  half  of  a  20  ns  rectangular  pulse.  While  a  drops 
closer  to  zero  after  10  ns  it  has  dropped  so  far  from  its  initial  value 
In  10  ns  that  very  little  further  €  is  generated.  While  examining 
the  data  in  Figs.  1  through  5note  that  the  assumed  temperature  is  quite 
high  so  that  a0=  0.508.  CARS  was  used  for  curves  A.l. 

Figure  2  shows  the  effect  of  keeping  the  CARS  term  only  in 
computing  A,  assuming,  when  calculating  x»  that  A_  and  AT  are  constant 

a  pi 

but  nevertheless  allowing  approximately  for  change  in  a  by  use  of 
Eq.  (20).  The  £  and  a  for  this  artificial  case  is  labeled 

a 

A. 2  and  is  computed  via  CARSO.  The  A  variation  is  only  approximate 
because  Ap  and  Aj  were  kept  constant.  Comparing  A. 2  to  a  correctly 

computed  curve  A. 3  and  to  the  45°  line  (for  CARS  only  and  a  =  A0)  gives 

an  idea  of  how  much  of  the  saturation  effect  is  due  to  a  changing  as 

opposed  to  how  much  is  due  to  the  neglect  of  the  stimulated  noncoherent 

generation  of  E  . 

a 

Figure  3  illustrates  the  effect  of  using  three  vibrational 

levels  instead  of  only  two;  The  idler  wave  Aj  in  the  three  level  case 

3 

is  taken  in  curves  A. 5  to  be  monochromatic  so  that  Ajj  *  10  , 

3 

Aj2  =  0.  For  curves  A. 4  Aj^  =  Aj2  =  10  . 

The  two  level  corresponding  curves  are  given  by  curve  A. 3. 

To  Interpret  the  curves  for  A. 5  a2  actually  becomes  negative  by  virtue 
of  stimulated  emissions  from  level  1.  We  have  plotted  |a2|  using  a 
dashed  curve  for  the  portion  of  the  curve  which  Is  plotted  when  a2 
Is  negative.  In  some  Instances  we  have  also  use<i  the  symbols ©and0 
to  indicate  the  sign  of  A;  or  a2. 
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1s  negative.  We  do  not  have  enough  data  to  draw  In  the  curve  In  the 
missing  Intervening  section  corresponding  to  the  I*Ij  range  In  which 
a2  goes  through  zero. 

We  note  that  although  the  energy  level  populations  are  quite 
radically  changed  by  the  more  complete  analysis  for  curves  A. 4  and  A. 5 
the  £  curves  of  A. 3  and  A. 5  differ  only  by  a  relatively  small  uniform 

A 

translation  for  all  I*IT  values.  The  £  plotted  for  curve  A. 4  is  the 

pi  3 

sum  of  £al  +  ^2*  It  corresponds  to  the  energy  a  receptor  of  uniform 
sensitivity  across  a  frequency  band  covering  n»al  and  u>a2  would  receive. 

Figure  4  presents  several  comparisons.  Curves  A. 3  vs.  A. 6  exhibit 
the  effect  of  Gaussian  pulses  (curve  A. 6)  as  compared  to  a  rectangular 
pulse  (A. 3).  The  two  pulses  have  the  same  duration  and  peak  power  so 
that  the  incident  Gaussian  pump  and  idler  pulses  deliver  81%  as  much 
energy  as  the  rectangular  pulses.  81%  of  the  <?  output  for  the  rectangular 
pulse  is  shown  as  an  additional  curve  in  the  figure. 

Curves  A. 3  vs.  A. 9  exhibit  the  effect  of  achieving  the  same 
range  of  l£lj  by  multiplying  Ij.  by  10  and  dividing  Ip  accordingly  by  /TO. 

The  large  difference  in  results  indicates  inaccuracy  in  the  simple  CARS  formula. 

In  Figure  5  the  curves  A. 3  vs.  A. 10  exhibit  the  effect  of 
changing  t  from  t  *  0.4  E  -  06  for  (A. 3)  to  t  *  0.1  E  -  06  for  A. 10, 
keeping  all  other  parameters  constant  (even  though  x  is  a  function  of, 
for  example,  temperature) .  The  difference  in  output  is  less  than  the 
accuracy  with  which  the  figures  are  drawn. 


Curves  A. 3  vs.  A. 7  show  the  effect  of  decreasing  the  coherence 
length  by  setting  Ak  *  Dk  *  6k  «  10.7  for  A. 7  against  setting 
Ak  =  Dk  =  5k  =  0  for  A. 3. 

Figure  6  depicts  the  direct  analog  of  Figure  3  but  carried  out 
for  the  data  of  Case  B  as  given  above  in  Section  C.  In  this  case  the 
excitation  A^  3  10  ,  Aj2  =  0  corresponds  to  curve  B.3  and  A^  =  Aj2  = 
to  curve  B.3.  The  two-level  case  is  curve  B.2. 


APPENDIX 


Boundary  Conditions 

To  solve  the  equations  (8}  and  (9)  requires  boundary  conditions 

on  A  and  F  *  (aA  /ax)  at  x  =  0. 
a  a  a 

For  the  Ao  we  have  Ap  and  Aj  given  nonzero  values  (determined  by 
the  experimental  set-up)  while  A.  and  Aw  *  0.  A  =  0  assumes  no 
anti-Stokes  is  generated  in  the  gas  outside  the  test  region— or  at 
least  that  It  is  negligible. 

To  obtain  boundary  conditions  on  the  derivatives  is  more 
difficult  in  a  sense  since  these  are  not  normally  specified  or  measured. 
We  use  the  following  approximate  procedure.  Referring  to  Eqs.  (8) 
through  (13),  we  first  consider  x's  independent  of  x  and  Fq  a  if^; 
i.e.. 


Then 


aA 


a  _ 


3X 


=  if 


3F 
_ a 

3X 


x=o 


a2A 


iaF 
_ a 

3X 


x=o 


x*o 


(A.l) 


(A. 2) 


where, in  evaluating  (3fa/3x)| x*0»we  use  the  boundary  conditions  on  A^ 

and  the  fact  that  F  ■  if  .  The  result  Is  that  the  ( 3 F ' /ax)Lso  can  be 
a  a  a  x~u 

written  explicitly  as  functions  of  Ap(0,t)  andA$(0,t).  Then 
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a 1  X=0 


.  J_L_  !fa  +  f  1 

[4,rv0  3*  “} 


(A. 3) 


These  are  the  values  used  to  start  the  solutions  of  (8)  and  (9). 
Specifically,  one  has 


jJ-  “  -12ir2vaxA2As  j  12n2C2x\)p[As[2  +  (v&x  +  v$x  )|Ap|2]  +  ak 


3^J  *  (12w2)2  v$AsI  ApJ 2  X*|va|Ap(2  x*  -  21|  As|  2vpXi  -  |Ap|2  v^x* 

0  L  j  1 

■  02it2)2  X Ap I As | 2  vp  j^jApI2  (1vsXi  -  vaxr)  -  vpX|As|2|  ^ 


where  x  *  xr  +  1x^  and  CSSRS  has  been  neglected. 
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